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We present an investigation of heat transport in gapless graphene-based Ferromagnetic /singlet Supercon- 
ductor/Ferromagnetic (FG|SG|FG) junctions. We find that unUke uniform increase of thermal conductance 
vs temperature, the thermal conductance exhibits intensive oscillatory behavior vs width of the sandwiched 
s-wave superconducting region between the two ferromagnetic layers. This oscillatory form is occurred by 
interference of the massless Dirac fermions in graphene. Also we find that the thermal conductance vs ex- 
change field h displays a minimal value at h/Ep ~ 1 within the low temperature regime where this finding 
demonstrates that propagating modes of the Dirac fermions in this value reach at their minimum numbers 
and verifies the previous results for electronic conductance. We find that for thin widths of superconducting 
region, the thermal conductance vs temperature shows linear increment i.e. To: T. At last we propose an 
experimental set-up to detect our predicted effects. 



I. INTRODUCTION 

In the recent years monoatomic graphite layers and 
structures containing this layers has attracted so much 
attentions theoretically and experimentally to itselP^^. 
The monoatomic graphite layer is called graphene and 
naturally treated as a two dimensional system. This 
t wo dim ensional artificial system made by Novoselov et 
(jJi3lTll_ jji graphene, low-excitation electrons follow the 
Dirac equation for their behaviors in various conditions 
and consequently Dirac equation can predict proper- 
ties of the systenJl^. Graphene exhibits very interesting 
properties that from one side confirm some of the pre- 
dicted phenomena in relativistic quantum mechanics e.g. 
Klein's paradox and from other side its high mobility and 
controllable Fermi energy in experiment make it very in- 
teresting and suitable in laboratory and industrj^illill^ 
From application point of view, it is very important to 
know the transport properties (charge, spin and thermal 
transport properties) of the devices including graphene 
junction&i^ 20--2J 

Thermal and charge transport properties are so much 
related to each other. Charge conductivity of the normal- 
superconductor (N/S) junctions at first was discussed by 
Blonder, Tinkham, and Klapwijk (BTK]P^ theoretically 
and their results could show good consistence with exper- 
iments. B-T-K take into account the contribution of the 
Andreev refiection'^ in the electronic transport of N/S 
junctions and use Boguliobov-de Genne (BdG) formal- 
ism to obtain the charge conductance at low tempera- 
tures. The BTK theory is limited to the clean regime of 
the heterojunctions, while cases with high impurities are 
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Heat Current Through the Junction 



FIG. 1. (Color online) Schematic model for a system of 
the sandwiched graphene-based s-wave superconducting layer 
with width L , between two uniform graphene-based ferro- 
magnetic substrates (FG/SG/FG). The ferromagnetic layers 
have similar directions in both sides. 



not withi n the regime of validity. Bardas et al. and Devy- 
atov et al^^^ generalized the BTK model for the charge 
transport and calculated the thermal current through the 
N/S junctions. For metals at low temperatures the ther- 
mal conductivity F has a linear behavior with respect to 
temperature i.e. F oc T, and their electronic conductivity 
reaches to constant value. So the Wiedemann-franz law 
is satisfied for thenPH. By depositing a superconducting 
electrode on a graphene substrate, the superconducting 
correlations can leak into the graphene, due to the prox- 
imity effect'^. Also ferromagnetism can be induced into 
the graphene by doping or using an external field^^ ^'^ 

In this paper we utilize the Dirac Bogolubov de-Genne 
equation and solve it for two dimensional systems includ- 
ing FG/SG and FG/SG/FG graphene junctions in which 
SG stands for the graphene-based s-wave superconduc- 
tor and FG stands for graphene-based ferromagnetic sub- 
strate. We use the generalized BTK formula for heat 
transport through the junctions and investigate the ther- 
mal transport properties of the FG/SG and FG/SG/FG 
junctions and as a especial case, we assume the exchange 
field h, equal to zero for approaching to normal case i.e. 
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NG/SG and NG/SG/NG. We find that the thermal con- 
ductance of the junction exhibits an intensive oscillatory 
shape, damping simultaneously by increasing the width 
of the superconducting layer. We find that the thermal 
conductance shows an exponential increase vs. temper- 
ature for large widths of superconducting layer that re- 
flects the s-wave symmetry of the Dirac ferniions inside 
the graphene as was mentioned by BTK but we find that 
for thin widths of superconducting layer within the tun- 
neling regime, thermal conductance T is linearly propor- 
tional to temperature T i.e. T ocT. Also we find that the 
thermal conductance vs. strength of the exchange field of 
the ferromagnetic substrate has a minimum near h Ep 
which this value moves toward smaller values h < Ep by 
increasing the temperature. 



II. THEORY 

We consider a graphene-based ferromag- 
netic/superconductor junction FG/SG which is placed in 
the a;?/— plane, and ideal interfaces between ferromagnet 
and superconductor located at a; = 0, are perpendicular 
to the X— axis. For investigating properties of the men- 
tioned system , one should solve the Di rac-B ogoliubov-de 
Gennes equation by following the Refs.'^^ 



Hq~ ah A 
A* -{Ho -ah) 



, (1) 



in which i?o (r) — —ihvp{arcdx+a-ydy)-\-U{r) — Ep, where 
Ox. and (Ty are 2x2 Pauli matrices. Also, A stands for 
order parameter of the superconducting layer, and one 
can mention the order parameter of the system under 
consideration as A(r, T) = A(T)0(x), in which Q{x) is 
the well known step function, Co- is the excitation en- 
ergy of the Dirac fcrmions with respect to the Fermi 
level, h(v) = hQQ{—x) is the exchange field energy of 
the ferromagnetic layer. Here cr = ±1 stands for spin-up 
and -down quasiparticles, and a — —a. U{r) shows the 
Fermi mismatch vector(FMV), and since through out the 
paper we investigate heavily doping cases, we assume a 
large value for mismatch potential in comparison with 
the Fermi energy Ep within the ferromagnetic region i.e. 
U{t) = -Uo6'(x) and Uq > Ep. Both of the and 
Vg- include two components of sublattices in hexagonal 
lattice of graphene. Thus each spinor in Eq.Q involves 
four components, and for the electron-like excitations in 
the ferromagnetic region (x < 0) we obtain: 



Y^COSQ^ 




, (2) 



and for hole-like quasiparticles: 
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(3) 



where fce(;j)_o- is component of the electron (hole)-like 
wave-vector, perpendicular to the interface and q is par- 
allel component of the wave vector which remains conser- 
vative during the scattering process. The above appeared 
factors 1 / ^/cos{a) and 1 / ■\/cos(a') guarantee same par- 
ticle current transport by the four wave-functions'^. Also 
Q!cr(a'g.) are injection angles of electron (hole)-like quasi- 
particles with respect to the axis normal to the interface 
(x-axis). They are defined as: 
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In the superconductor region (x > 0), wave function for 
the hole-like quasiparticles read as: 
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and for the electron-like quasiparticles: 



(8) 



± 

S,e 



J{±iko+ix)x+qy) 



1 



where 



-icosh^i(^),e < Ac 



Uo + Ep 



Uo 



hvp 

Ep 



ko{hvpy 



sin /3, 
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(11) 
(12) 



(13) 



Uo + Ep' 

here vp is energy-independent Fermi velocity in 
graphene. We define right going (-|-a;-direction) and left 
going (— x-direction) quasi-particle wave-functions with 
plus and minus signs e.g. ip~^ , ip~ respectively. In the 
mean field approximation, we assume high doping regime 
Uo+Ep » AcP. It is clear that during the scattering pro- 
cess from the interface, q component of parallel wave vec- 
tor and the energy of quasiparticles are constant (elastic 
scattering). The wave functions of the moving quasiparti- 
cles must satisfy the boundary conditions at the interface 
between ferromagnet and superconductor as in Ref.'^. For 
FG/SG junction the boundary condition reads as: 
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FIG. 2. (Color online) The normalized thermal conductance 
of single FG/SG junction vs temperature for three values of 
the ferromagnetic sheet's exchange field. 



FIG. 3. (Color online) The normalized thermal conductance 
of the single FG/SG junction vs exchange field of the ferro- 
magnetic layer for three values of temperatures . 



For the FG/SG/FG structure that schematically is shown 
in Fig. [T] however the boundary conditions at a; = and 
X — L respectively are: 



(15) 



and 



(16) 

Here, r^.o- is amplitude of the Andreev reflection, rjv.o- is 
amplitude of the normal reflection, t^^a and tfi,a are am- 
pHtudes of the electron-like and hole-like quasiparticle's 
transmission, respectivley. Substituting the wave func- 
tions into the above boundary conditions, we calculate 



the coefficients in Eq. (jM]), for FG/SG and FG/SG/FG 
junctions. Then we now are able to obtain the probabil- 
ity of the Andreev reflection {RA,a = k^.o-P) and normal 
reflection {Rn,<7 — ^A^.trP)- As seen in Fig. [l] interfaces 
are normal to the x-axis and superconductor region is 
between a; = and a: = L, so the x-dcpcndent order pa- 
rameter can be written as A(a;) = Af)d{x)6{L — x). 

The nor maliz ed thermal conductance T = F /Fq is given 
as follo^^^2SI3Il 



r' ^Fo^ / / dEdaaCos{aa){l- \ rN,a{E,aa) \' 



T2 C0sh2(J.)' 



(17) 



where Fq = Ep /27T^h?VFkBAo is a constant. Here all 
parameters are normalized, i.e. energies, with respect to 
pair potential at zero temperature (Aq = A(r = 0)) , 
and temperatures, with respect to the critical tempera- 
ture of superconducting order parameter (Tc) . Through- 
out the paper we set Ag = h — kg ^ ^ ui our computa- 
tions. 



III. RESULTS AND DISCUSSION 



Now we proceed to present the main results of the 
paper. In order to obtain more realistic experimental 
results, we consider highly doping junctions and set a 
large mismatch potential Uq. The Fermi energy of the 
graphene is externally controllable and can be tunecP'^. 
Energy of quasiparticles remains within our regime of 
validity, for values near to leV i.e. the quasiparti- 
cles follow the Dirac equation. Throughout the paper we 
fix the Fermi energy Ep = IOAq which typically places 
Fermi energy within 10-15 meV. We change the width of 
the junction from L ~ ^5 up to L ~ 9fs and investigate 
how the thermal conductance is varied by the variation 
in L/S^s- Also we investigate other possibilities of varia- 
tions in every parameter involved the problem and how 
they infiuence the thermal conductance. When we need 
fixed temperature, use T /Tc = 0.2, when we need fixed 
exchange field use h/Ao = 8 and when we need a fixed 
width, typically use L/^s — 4. As will be discussed in 
the following sections, we find that unlike the exponen- 
tial increase of the normalized thermal conductance with 
respect to the temperature in the FG/SG junctions re- 
flecting the s-wave superconducting correlation, the nor- 
malized thermal conductance for FG/SG/FG junctions 
shows a linear increasing with respect to the tempera- 
ture for small widths (L ~ ^5) of the superconducting 
region i.e. F cx T. For FG/SG/FG junctions, the ther- 
mal conductance vs. width of the superconducting region 
L/^s shows an intensive oscillatory behavior, also it has 
a minimum vs. strength of the exchange field in both 
semi-infinite ferromagnetic graphene sheets. 



FIG. 4. (Color online) The normalized thermal conductance 
r vs the normalized width -L/^s of graphene-based supercon- 
ducting layer for different values of h in FG/SG/FG junction. 
The temperature is fixed at T/Tc = 0.2. 



A. Heat transport through the single 
Ferromagnetic/Superconductor junction 

At first, we investigate how the normalized thermal 
conductance of the single junctions (FG/SG) behavior 
by changing both the temperature and exchange field 
h of the semi-infinite ferromagnetic sheet. We fix the 
Fermi levels of superconductor and ferromagnetic sheet 
at Ep — 10 Aq and use a large mismatch potential t/o, 
thus with such parameters we remain in the heavily 
doped regime. As seen in Fig. 2, the thermal conduc- 
tance increases exponentially by increasing the temper- 
ature T, that verify the presence of the induced s-wave 
correlation in the graphene sheet. As seen in Fig [2] the 
curve of the thermal conductance for /i/Ag = 15 is placed 
between that of the thermal conductance for /i/Aq = 
and 9. This findings show that the thermal conductance 
has a minimum vs. h/Ao as seen in Fig. [sj In Fig. 
[3] the normalized thermal conductance reach to a mini- 
mum near /i ~ i^j?. Therefore by increasing the exchange 
field of the ferromagnetic layer, propagating modes of the 
Dirac Fermions decay , and reach to a minimum num- 
ber near Fermi level. By increasing the temperature, 
the order parameter of the superconducting correlations 
decays, thus increasing of the temperature helps the ex- 
change field to make propagating modes reach to their 
minimum value at smaller values of exchange field h as 
seen in Fig. [3j Now we proceed to present our findings 
for double FG/SG/FG junctions. 



B. Heat transport through the double 
Ferromagnetic/Superconductor/Ferromagnetic junction 

Now we turn our attention to the second graphene- 
based structure, namely double FG/SG/FG junction. 
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FIG. 5. (Color online) The normalized thermal conductance 
r vs temperature for three values of normalized widthes of the 
superconducting layer L/^ — 1, 4, 15 in FG/SG/FG junction. 
The exchange field is fixed at /i/Ao=8. 



The suggested set-up of the junction schematically is 
shown in Fig. [T] Such junctions can be made by deposit- 
ing a superconducting electrode on top of the graphene 
sheet. As before we fix the Fermi energy of the supercon- 
ductor and ferromagnetic at Ep = IOAq. As shown in 
Fig|4j the normalized thermal conductance shows inten- 
sive oscillatory behavior versus width of superconductor 
region. This manner is a consequence of the coherent in- 
terference of the Dirac fermions in facing with two inter- 
faces. Also, by increasing the width of superconducting 
region, the normalized thermal conductance loses the am- 
plitude of its oscillation and shows exponentially decre- 
ment towards single junction. When the Dirac fermions 
lose their coherency , they show the mentioned manner 
and this illustrates the fact that when the width of the su- 
perconducting region becomes larger, the probability of 
tunneling of the Dirac fermions across the superconduct- 
ing region is reduced. The thermal conductance versus 
temperature is shown in[5j The thermal conductance for 
small widths of the superconducting region shows a linear 
behavior similar to the behavior of the thermal conduc- 
tance of metals, discussed in the introduction. This man- 
ner reflects the tunneling process of the Dirac fermions 
through the superconducting region. Thus one can ex- 
pect that when the width of this region goes to larger val- 
ues, this phenomena (tunneling) is reduced and the ther- 
mal conductance behaves like the single junction which 
is clearly shown in Figj5] for L/^s — 15. The thermal 
conductance as a function of the exchange field of the 
ferromagnetic region is plotted in Fig. |6] The behaviors 
of the thermal conductance vs. h for double junction 
(FG/SG/FG) and the single junction (FG/SG) are the 
same. The thermal conductance reach to a minimum 
near h ~ _E^,and the minimum value is moved to smaller 
ones by the temperature increment. 
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FIG. 6. (Color online) The normalized thermal conductance 
r vs strength of normalized exchange field h/Ao of the ferro- 
magnetic sides for three values of temperatures in FG/SG/FG 
junction. The width of the superconducting layer is fixed at 
i/Cs=4. 



IV. SUMMARY 



In the present paper we have considered two 
types of gapless graphene-based junctions, FG/SG and 
FG/SG/FG with s-wave superconducting electrodes, and 
we have investigated the heat transport properties of the 
systems. In particular, we have investigated how the 
variations of the variable quantities in the system can 
influence the thermal conductance of the junctions. Uti- 
lizing the Dirac-Boguliobov de Genne equation for quasi- 
particles inside the graphene and appropriate boundary 
conditions for the obtained wave functions, we have de- 
rived the Andreev and normal reflection cocfflcicnts, and 
used them for calculating the normalized thermal con- 
ductance numerically. We found that for single junction 
(F/S) the normalized thermal conductance indicates the 
previous exponential form vs. temperature. Increasing 
exchange field h lowers the propagating modes of the 
Dirac Fermions down to a minimum value. Due to the 
tunneling phenomena of the Dirac fermions through the 
superconducting layer, the heat transport properties of 
the double F/S/F junctions are different from single F/S 
junction. The thermal conductance vs temperature ex- 
hibits linear behavior in the tunneling limit (small widths 
of the superconducting region). The thermal conduc- 
tance vs. width of the superconducting region exhibits 
very intensive oscillatory behavior and also the quantity 
indicates previous behavior of the single junction (F/S) 
vs. h i.e. the enhancement of the exchange field from 
up to a value near Ep, lowers the propagating modes of 
the Dirac fermions. 
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